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ABSTRACT 


Chandra X-ray imaging spectroscopy of the nearby Herbig Ae star HD 163296 at 100 AU angular 
resolution is reported. A point-like, soft (kT~ 0.5 keV), emission-line source is detected at the location 
of the star with an X-ray luminosity of 4 x 10 29 ergs s _1 (logLx/Tboi = —5.48). In addition, faint 
emission along the direction of a previously-detected Lya-emitting jet and Herbig-Haro outflow may be 
present. The relatively low luminosity, lack of a hard spectral component, and absence of strong X-ray 
variability in HD 163296 can be explained as originating from optically-thin shock-heated gas accreting 
onto the stellar surface along magnetic field lines. This would require a (dipole) magnetic held strength 
at the surface of HD 163296 of at least ~100 G and perhaps as high as several kG. HD 163296 joins 
the T Tauri star TW Hya in being the only examples known to date of pre-main-sequence stars whose 
quiescent X-ray emission appears to be completely dominated by accretion. 

Subject headings: stars: individual (HD 163296) — stars: emission-line, Be — stars: pre-main sequence 
— X-rays: stars 


1. INTRODUCTION 

The Herbig Ae/Be stars (Herbig 1960) are the more 
massive counterparts of T Tauri stars. They are character¬ 
ized by a strong IR excess, emission lines, and luminosity 
class III to V (Waters & Waelkens 1998). Herbig Ae/Be 
stars represent an important link between high mass stars 
that evolve directly from embedded protostars to the zero- 
age main sequence (the onset of core H-burning) and low 
mass stars that stop accreting before reaching the main 
sequence and that can be observed in this hydrostatic pre- 
main sequence phase (Appenzeller 1994). 

Herbig Ae stars are often strong X-ray emitters with lu¬ 
minosities of order 10 29 to 10 32 ergs s -1 (e.g., Hamaguchi 
et al. 2002; 2004). This high level of activity is possibly 
linked to the presence of an accretion disk and is probably 
ultimately driven by magnetic fields. However, the origin 
of magnetic fields in Herbig Ae stars is unknown. Main 
sequence A stars lack the thick surface convection zones 
(Gilliland 1986) of late-type stars (and their T Tauri pro¬ 
genitors) that drive magnetic dynamo activity and power 
X-ray-emitting chromospheres and hot coronae. Hence, 
they are weak X-ray emitters (Simon, Drake, & Kim 1995; 
Linsky 2003) or even X-ray dark (Schmitt 1997). Per¬ 
haps held generation in Herbig Ae stars is not through the 
solar-type a-tt dynamo mechanism but from compression 
of interstellar cloud fields in the initial collapse (Dudorov 
et al. 1989), induced by deuterium-shell burning during 
contraction towards the main sequence (Palla & Stahler 
1993), or generated via internal differential rotation dur¬ 
ing an early accretion phase (Tout & Pringle 1995). Such 
fields may decay rapidly through turbulent magnetic dif- 


fusivity as the radiative core develops and only be present 
in the youngest Herbig Ae stars. 

In analogy to the T Tauri stars, whose X-ray activity 
is also often much greater than in evolved stars of similar 
mass, perhaps circumstellar material allows other mag¬ 
netic configurations such as star-disk and disk-disk fields 
to generate the strong flares, hard X-ray emission, and 
high X-ray luminosity seen in some intermediate-mass Her¬ 
big Ae systems. Energetic reconnection flares may arise at 
the corotation interface between dipole stellar fields and 
the inner disk (Hayashi, Shibata, & Matsumoto 1996; Shu 
et al. 1994; Birk et al. 2000) or through differential rota¬ 
tion of the field-threaded disks themselves (e.g., Romanova 
et al. 1998). Magnetic fields are believed to regulate pro- 
tostellar collapse through ambipolar diffusion, to funnel 
accretion from the disk, and to transfer disk orbital mo¬ 
tion to collimated outflows (see Feigelson & Montmerle 
1999 for a review). 

HD 163296 is a nearby (122 pc) Herbig Ae star (AlVe, 
Au=0. m 25) with an effective temperature T e //=9300 K, 
luminosity log(L*/LQ)=l. 481 ^ 10 : mass Af* = 2.3M 0 , 
and radius I?* = 2.1Rq (van den Ancker, de Winter, & 
Tijn A Djie 1998). Based on its position in the HR dia¬ 
gram, HD 163296 is at an age t= 4+® 5 Myr (van den An¬ 
cker et al. 1998); intermediate between heavily-embedded 
pre-main-sequence A stars and near-ZAMS stars with de¬ 
bris disks. HD 163296 has a strong infrared excess (Hillen- 
brand et al. 1992; Meeus et al. 2001) and variable Balmer 
series emission lines (Baade & Stahl 1989; Pogodin 1994; 
Beskrovnaya et al. 1998) typical of the Herbig Ae class 
though it has no associated molecular or dark cloud (e.g., 


1 Universities Space Research Association, NASA Marshall Space Flight Center, SD50, Huntsville, AL, USA 

2 Smithsonian Astrophysical Observatory, MS-3, 60 Garden Street, Cambridge, MA, USA 

3 Space Science Department, NASA Marshall Space Flight Center, SD50, Huntsville, AL, USA 

4 Eureka Scientific; and Laboratory for Astronomy and Solar Physics, Code 681, NASA Goddard Space Flight Center, Greenbelt, MD, USA 

5 Thomas Aquinas College, Santa Paula, CA; and Institute for Astrophysics & Computational Sciences, The Catholic Univ. of America, 
Washington, DC, USA 

6 NASA Goddard Space Flight Center, Code 681, Greenbelt, MD, USA 

1 


2 


The et al. 1985). The large infrared excess arises from 
heated, optically thick dust within a circumstellar disk 
(Hillenbrand et al. 1992; Meeus et al. 2001). The disk 
has a radius of ~450 AU (3."7) viewed at ~60° inclina¬ 
tion (Mannings & Sargent 1997; Grady et al. 2000). The 
time-variable double-peaked and P Cygni Balmer and UV 
emission line profiles arise from a stellar wind, an extended 
chromosphere, and/or rotation (e.g., Catala et al. 1989). 
Radio continuum observations (Brown, Perez, & Yusef- 
Zadeh 1993) of HD 163296 also suggest wind-driven mass 
loss. Orthogonal to the disk is an axially-aligned chain 
of Herbig-Haro nebulae extending several seconds of arc 
above and below the disk (Grady et al. 2000). A Lya- 
bright jet extending 6" along the SW arm of the Herbig- 
Haro flow has also been discovered (Devine et al. 2000). 

HD 163296 was examined with Chandra in order to 
study the spatial distribution of the X-radiation from the 
star and circumstellar environment (§03 and to character¬ 
ize the X-ray spectrum and light curve (§ 0 . The com¬ 
bination of an accretion disk and jet in a relatively old 
Herbig Ae object, a low extinction due to a lack of as¬ 
sociation with dark clouds, the presence of Herbig-Haro 
objects along the jet axis, and an otherwise source-free 
field make HD 163296 a particularly promising object for 
high-resolution X-ray imaging spectroscopy. It is found 
that the bulk of the X-ray emission is associated with the 
central star and its immediate surroundings. This emis¬ 
sion can be described as a cool, ~0.5 keV, emission-line 
plasma with a luminosity ~4 x 10 29 ergs s -1 in the 0.3- 
3.0 keV range. There is also evidence for X-ray emission 
along the jet. Possible origins of the X-ray emission are 
discussed in § 4. 

2. X-RAY IMAGE 

HD 163296 was observed for 19.2 ks on 10 August 
2003. The observation utilized Chandra’s Advanced CCD 
Imaging Spectrometer (ACIS) operating in 1/4-subarray, 
timed-exposure, imaging mode. Only the central three 
CCDs, S2-S4, were active during the observation result¬ 
ing in an approximately 2/1x25/3 image and a 0.941 s 
frame time. This configuration was chosen to provide the 
highest integrity image by minimizing event pileup. Ran¬ 
domization of event locations at the sub-pixel level, part 
of the standard data processing pipeline, was removed by 
reprocessing the data. Standard ASCA grade 02346 events 
in the 0.3-8.0 keV energy range were included in the anal¬ 
ysis. 

Nine weak point-like sources (each with 9 to 17 counts 
detected corresponding to signal-to-noise ratios of 2.5 to 
3.6) and one bright source (1131 counts) were detected 
in the held. One of the weak sources is 0."20 offset in 
Right Ascension from a source in the Two Micron All 
Sky Survey (2MASS) All-Sky Data Release catalogue. 
This is comparable to the estimated statistical uncertainty 
in the X-ray position of that source. No other astro- 
metric catalog matches were found for any of the other 
weak X-ray sources. The bright source is offset -0."26 
in Dec, 0.00 s in RA, from the Hipparcos position of 
HD 163296 extrapolated to the time of the X-ray ob¬ 
servation (17 h 56 m 21/28524, -21 d 57 m 22/0248). This 
is within the la nominal Chandra absolute pointing un- 

7 http://cxc.harvard.edu/cal/ASPECT/celmon/ 


certainty 7 of 0."38 and infers the X-ray emission is from 
HD 163296 and not from a nearby late-type companion. 

Hubble observations also place tight constraints on a 
potential companion. Late-type companions produce con¬ 
spicuous diffraction patterns in STIS chronographic im¬ 
agery with spatial intensity profiles distinctly different 
from those of early A stars (see Grady et al. 2004a for 
examples in the Held of HD 104237). STIS imagery of 
HD 163296 shows no evidence for a late-type companion 
beyond 0."05 of HD 163296. In addition, long-slit STIS 
spectra of HD 163296 are characteristic of chromospheric 
and transition region spectra of early-type stars and dis¬ 
play the spatial profile of an unresolved point source. The 
emission line features seen in these and in FUSE data 
(Grady et al. 2004b) are uniformly broad and often double- 
peaked. Such features are characteristic of Herbig Ae stars 
and are much broader than those of late-type stars; even 
rapid rotators like the MlVe star AU Mic (Pagano et al. 
2000; Robinson et al. 2001). Both line and UV/FUV emis¬ 
sion are spatially coincident, at the resolution of Hubble , 
with the photosphere of HD 163296. In particular, emis¬ 
sion associated with transition region temperature plasma, 
N V and O VI, originate on the rapidly rotating HD 163296 
rather than on a more slowly rotating, but otherwise un¬ 
seen, late-type companion. We conclude there are no late- 
type stellar companions to HD 163296. 

The central 40" x 40" region of the X-ray image is shown 
in the left panel of Figure 1. A nearest-neighbor (boxcar) 
smoothing has been applied to the image and contours 
have been added to emphasize the weak extension from the 
source toward the SW and diffuse structure located ~7" 
from the bright source along this jet axis. The roll angle 
for the observation is 268.°5 placing any readout streak 
along an east-west line (left-right in the figures) through 
the bright source or roughly 45° from the jet axis. Thus, 
there is no contribution to the X-ray extension along the 
jet axis from the readout streak. The center panel of Fig¬ 
ure 1 shows the positions of the individual X-ray events 
with the location of the weak diffuse feature marked. 

The right panel of Figure 1 shows, for comparison, a 
portion of an ultra-narrowband [SII] image of HD 163296 
taken on 2004 May 17 using the Goddard Fabry-Perot 
(GFP) instrument and a coronagraphic wedge to occult 
the central starlight on the 3.5m telescope at the Apache 
Point Observatory in Sunspot, New Mexico. The image 
was constructed from an on-band image centered at 6724 
A with a 12 A FWHM ( or 540 km/s, thereby sampling 
both the jet at -300 km s -1 and the counterjet at +340 
km s -1 ) from which an off-band image centered at 6750 A 
with 12 A FWHM was subtracted. Both images are 900 s 
exposures taken under conditions of l."5 seeing. The GFP 
Held of view is a circular aperture of 3/7 diameter and the 
CCD plate scale is 0."366 pixel -1 . A 40"x40" subim¬ 
age of the full GFP image centered on the star is shown. 
Herbig-Haro knots associated with the jet and counter jet 
of HD 163296 are visible in the raw GFP observation, and 
can be traced out to ±27" from the star after subtraction 
of the oH-band image. The knots in both the jet and coun¬ 
terjet were measured at S/N=10. The weak X-ray feature 
does not overlap any of the optical knots when their po¬ 
sitions are devolved to the time of the X-ray observation 
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Fig. 1.— Left panel: Smoothed 40"x40" ACIS image of the region around HD 163296 in the 0.3 to 3.0 keV bandpass obtained in 2003 
August. Tick marks denote steps of five ACIS pixels (0. ,, 492 pixel -1 ). North is up and East is to the left. The star symbol marks the 
Hipparcos position of HD 163296. Countours denote 0.1, 0.2, 0.5, 2, 5, 10, and 30 c pixel -1 levels. The lowest contour extends >2" from 
the center of the bright source along the SW in the direction of the Ly a jet reported by Devine et al. (2000) (P.A. 42.°5±3.°5). The diffuse 
isolated feature is ~7" from the central source and nearly along this axis. The debris disk (Mannings & Sargent 1997; Grady et al. 2000) is 
oriented perpendicular to the jet with radius ~450 AU (3."7) and inclination ~ 60°. Center panel: Individual X-ray event locations for the 
data shown at left. The symbol denotes the centroid of the weak diffuse feature. Right panel: 20" x20" continuum-subtracted [S II] image 
obtained 2004 May. Note the occulting wedge shadow and diffraction spikes from HD 163296. The jet and knot structures are clearly visible 
and extend ~27" from the star along both the jet and counterjet. The structure in the counterjet (toward the NE) has no apparent X-ray 
counterparts. The optical knots have proper motions of 0.410±0.005 sec yr -1 . 


using the known knot proper motions. 

The spatial distribution of the X-ray events from the 
central portion of HD 163296 was compared to a model 
point spread function (PSF) to test for possible source ex¬ 
tension. The model PSF was generated using the Chandra 
Ray Tracer (ChaRT) web-based interface 8 to the SAOsac 
raytrace code using the best-fit singl e-kT model spectrum 
(kT ~ 0.5 keV, § 3). The radial profile (averaged over az¬ 
imuth) of the data shows a slight excess above the model 
PSF in an annulus between 0.5 and ~2."5 of the emission 
peak. However, the ray-trace software does not account 
for spacecraft dither and the observed excess is consistent 
with this model limitation. Therefore, except for emis¬ 
sion along the jet, there is no conclusive evidence that the 
bright source is extended. 

3. X-RAY SPECTRA & TIMING 
3.1. HD 163296 

A time-dependent gain 9 correction was applied to the 
event file. The spectrum was then extracted from a 2- 
pixel-radius region about the center of the bright source. 
Corresponding redistribution matrix and auxilliary re¬ 
sponse files were generated using CXC/CIAO tools (v.2.3). 
Models were fit to the spectrum using the XSPEC (v.11.2) 
spectral-fitting package. The acisabs multiplicative ab¬ 
sorption model 10 was included to account for the temporal 
decrease in low-energy sensitivity of ACIS. The C-statistic 
fit-minimization algorithm was used to determine the best- 
fit parameters for the unbinned data. 

Figure 2 shows the spectrum, model, and fit residuals for 
an absorbed variable-abundance (vmekal) emission-line 
model (C-statistic 214.5 for 184 spectral bins). For display 
purposes, up to 5 channels have been combined to achieve 

8 http://asc. harvard.edu/chart 

Al I p://hea-www.harvard.edu/alexey/acis/tgain 

1 °ht Ip://www. astro.psu.edu/users/chartas/xcontdir/xcont.html 


a minimum 5cr significance. The parameter estimates 
and 90% confidence ranges for this model are: hydro¬ 
gen column density, A^///10 20 =7.62 ± 1.85 cm -2 ; plasma 
temperature, fcT=0.49 ± 0.03 keV; and emission measure, 
EM=(1.1 ±0.3) x 10 53 cm -3 . The metal abundances were 
allowed to vary in the model fits but were restricted to pre¬ 
serve solar ratios for certain element groups whose first ion¬ 
ization potentials (FIPs) are similar (to allow for possible 
FlP-based chemical fractionation commonly observed in 
late-type stellar coronae) and whose principal X-ray lines 
are formed at similar energies; namely Fe-Ni-Ca, Ne-Ar, 
Mg-Si-S-Al-Na, and C-N-O. In this way, the number of 
free parameters is not excessive. The estimated best fit 
metal abundances relative to solar are 0.13±0.02 for the 
Fe group, 0.05±0.05 for the Ne group, 0.18±0.08 for the 
Mg group, and 0.25±0.12 for the C group. 



0.5 1 2 

channel energy (keV) 











4 


Fig. 2. — Upper panel: Observed spectrum of HD 163296 (sym¬ 
bols) and best-fit single-temperature thermal emission-line (vmekal) 
model (solid line). Lower panel: Fit residuals. 

The best-fit value of Nh is slightly higher than esti¬ 
mated from the optical extinction, Ay= 0.25 mag, assum¬ 
ing a normal extinction law (e.g., Bohlin, Savage, & Drake 
1978). The model luminosity is 2.9 x 10 29 ergs s -1 corre¬ 
sponding to an intrinsic (absorption-corrected) luminosity 
of L x =4.0 x 10 29 ergs s -1 in the 0.3-3.0 keV energy range 
after correcting for the finite spectral-extraction region. 

Additional models were applied to try to improve the 
fits and to better characterize the composition and intrin¬ 
sic absorption column. These included relaxing the metal 
abundance restrictions in the original model, 2- or 3- 
temperature vmekal models, and variable-abundance ab¬ 
sorbing columns. None of these models gave statistically 
significant improvements. The temperature and overall 
abundance are constrained by Fe which emits numerous 
strong Fe XVII lines in the ~0.7 to 1.0 keV band. A lower- 
temperature plasma (in ionization equilibrium) would be 
needed to produce a stronger O VII triplet at ~0.57 keV 
and H-like N Lya at 0.5 keV; merely increasing O and 
N abundances does not alter the model spectrum. If the 
abundances of Mg and Si are ungrouped and increased to 
~0.2 Zq, then the emission at 1.35 keV (He-like Mg) and 
1.84 keV (He-like Si) increases but the fit is not statisti¬ 
cally improved. Ne X Lya also contributes at 1 keV but 
the best-fit Ne abundance is slightly lower than for other 
elements at 0.05 times the solar photospheric value. 

We also investigated models with an additional thermal 
component, though these models resulted in no significant 
improvement to the model fit for the case of both lower 
and higher temperatures. Forcing a component at a fixed 
temperature of 3 keV resulted in a contribution of at most 
10 % to the total flux. 

Accretion shocks or shocks from a stellar wind might 
contribute to the observed X-ray emission from HD 163296 
(§ EUJ- Models based on the physics of shock-heated 
plasma (e.g. pshock and variants following the work of 
Borkowski et al. 2001), were also fit to the spectrum. The 
resulting fits were no better than the simple absorbed 
vmekal model. The best-fit parameters of the shocked 
gas models (absorbing column, post-shock temperature, 
abundance, and emission measure) were comparable to the 
values determined from the emission-line model. 

3.2. X-ray Knot 

There are 5 counts in the 0.3-3.0 keV band detected 
within a 2."5 radius around the position of the “knot” 
shown in Figure 1. This feature was not detected with 
our source-finding algorithm. However, the expected back¬ 
ground is only 0.6 c giving a probablity of 0.0004 that 5 
or more Poisson-distributed counts would be detected in 
a given region of this size. On the entire S3 chip (in 1/4- 
subarray mode) there are 7 of 2625 regions of this size 
found to contain 5 or more counts. Thus the distribution is 
not strictly Poissonian; there are also real sources present 
below the detect limit. There are 13.75 regions of this size 
within the area of the optical jet structure extending ±27" 
from HD 163296. Within this area, the probability of de¬ 
tecting a 5-count feature in a 2."5 radius circle is therefore 
0.04. 


The source counts are clustered near 0.8 keV, corre¬ 
sponding to the peak of the spectrum of a thermal plasma 
of temperature ~0.3 keV when folded with the Chandra 
response. The corresponding unabsorbed X-ray luminos¬ 
ity of the knot is ~ 3 x 10 2 ' ergs s -1 . 

3.3. HD 163296 Light Curve 

The 20 ks Chandra observation corresponds to about 
1/3 of the rotation period of HD 163296 (van den Ancker 
et al. 1998). No flaring of the source was observed over the 
short observation. The binned light curve of HD 163296 
(Figure 3) is consistent with a constant flux (y 2 =31.9 for 
25 dof) with variations not more than ~20% in amplitude. 
A Kolomogorov-Smirnov statistic was also computed to 
test the source for time variability by comparing the cu¬ 
mulative event arrival times binned at the frame time of 
0.94 s to that expected for a steady source. This test 
showed the source to be marginally variable (significance 
level 0.006) because of the shallow dip in the light curve 
about mid-way through the observation. These two tests 
indicate there is no evidence for strong variability during 
the observation. 



Fig. 3. — X-ray light curve of HD 163296 accumulated into 800 s 
bins. Data was extracted from a 2-pixel radius about the peak of 
the spatial distribution of events. Only events in the 0.3—3.0 keV 
energy range are included. The background contributes less than 
one event to the >1000 events in the light curve. 

4. DISCUSSION 
4.1. HD 163296 

The Chandra observation reveals a point-like object 
within 0."25 (30 AU) of the Hipparcos position of 
HD 163296. The X-ray luminosity is 4xl0 29 ergs s" 1 
(logTx/Tboi = —5.48) and the best-fit temperature is 
0.5 keV (6xl0 6 K). Notably, there is no evidence for 
highly-variable, hot (fcT>3 keV), and luminous (up to 
~10 32 ergs s -1 ) X-ray emission that has been attributed 
to reconnection flares in other Herbig Ae systems (e.g., 
Hamaguchi et al. 2002; 2004) in analogy to the late- 
type T Tauri stars. The spatial coincidence, X-ray spec¬ 
trum, and lack of X-ray flares strongly suggest the X-ray 
emission is associated with HD 163296 and not with a 
hidden late-type companion. This conclusion is further 
supported by the lack of spatial and spectral evidence 
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from Hubble /STIS imagery for a late-type companion to 
HD 163296. 

Line-driven winds also cannot account for the 
X-ray emission. The wind mass-loss rate is 

M~lO~ 15 (L*/L 0 ) 2 (Af*/M Q ) -1 M q yr -1 (see Abbott 
1982). For HD 163296, 5.5 x 10 “ 13 M 0 yr " 1 

and there is too little kinetic power in the wind, 
£’=(1/2 )Mu 2 ~6 x 10 28 ergs s -1 , to account for the X- 
ray luminosity even if all this energy were converted to 
X-radiation through shocks. 

The presence of a circumstellar debris disk, jets, and 
Herbig-Haro outflows are all strong indications of accre¬ 
tion and accretion-induced outflow in HD 163296. The 
luminosity released in accretion is L acc =GM+M/2R ir (ex¬ 
cluding the portion viscously dissipated in the disk). For 
HD 163296, L a cc =6.6xl0 33 A'f 7 ergs s _1 where AL 7 is the 
mass accretion rate in units of 10“' M 0 yr -1 . 

In the absence of strong magnetic fields, the accretion 
energy is released in an optically-thick boundary layer at 
the surface of the star. In this case, the accretion en¬ 
ergy radiates as a blackbody from an area ~ 47 ri?*L?, where 
IL~0.1I?* is the height of the emission region (Popham et 
al. 1993), with a temperature T^{GM i ,M/SuRZ/Ha) 1 ^ or 
~7900K for HD 163296. Boundary layers will there¬ 
fore be copious UV emitters but cannot account for the 
X-ray spectrum of HD 163296. 

Alternatively, the inner regions of the disk could be 
disrupted by a magnetic field associated with the accret¬ 
ing star. In this case, disk material couples to the field 
lines and is channeled onto small regions near the mag¬ 
netic poles at an angle nearly normal to the stellar surface 
(Konigl 1991, see, e.g., von Rekowski & Brandenburg 2004 
for a more recent study). Accreted material is strongly 
shocked on impact. Applying the jump conditions for a 
shock perpendicular to the flow implies a post-shock tem¬ 
perature T s = 2>GM i ,pm p /&kR+ where n=0.6 for ionized 
solar abundances, m p is the proton mass, and the velocity 
just upstream of the shock is approximately the free-fall 
value, n//=(2GAL*/L?*) 1 / 2 , at the surface. For HD 163296, 
T s ~5.7x 10 6 K — consistent with the observed X-ray tem¬ 
perature. 

If the accretion rate is high, however, then the stream 
can be of sufficiently high density that it will penetrate 
the stellar photosphere. Much of the X-radiation is then 
attenuated by preshock gas within the accretion flow or 
reprocessed within the stellar photosphere (e.g., Ulrich 
1976; Stahler, Shu, & Taam 1980; Lamzin 1995; Calvet & 
Gullbring 1998; Lamzin 1999; Drake 2004). The radiation 
field will become optically thick with an effective temper¬ 
ature of T~(GA/*M/ 47 ri? 3 /cr ) 1 / 4 where / is the fraction 
of the surface area covered by the accretion stream. For 
/ = 0.06 (Shang et al. 2002), the effective temperature for 

HD 163296 is only ~6000M 7 1/4 K. 

The depth in the atmosphere at which the accretion 
shock will form is given by the balance between the ram 
pressure of the stream, p ram = 2 Afv/// 47 ri? 2 /, and the 
photospheric gas pressure. For typical gas pressures of a 
few 10 2 - 10 3 dyne cm -2 , the mass accretion rate must be 
below 0.1 for the X-ray emitting shock to lie either 
above or only partially within photospheric layers. For 
significantly lower filling factors of less than 1 % favoured 


by some analyses of accreting T Tauri stars (e.g. Calvet 
& Gulbring 1998), the X-ray shock can remain completely 
buried in the photosphere even for accretion rates as low as 
M 7 ~0.01. Thus, the observed X-rays from HD 163296 can 
only be consistent with an accretion scenario if the shocked 
material lies above, or at very shallow depths within, the 
photosphere and this occurs only for accretion rates below 
~10 “ 8 M 0 yr” 1 . 

The lower limit to the accretion rate obtained by assum¬ 
ing all the accretion luminosity emerges as X-radiation is 
only M 7 ~6x10 “ 5 or lif~6xl0 “ 12 M 0 yr -1 . At such a low 
rate most of the accretion luminosity should be observed 
in X-rays for plausible ranges of filling factor. The actual 
accretion rate, however, must be higher as only a fraction 
of the available energy emerges in X-rays and a portion of 
the accretion energy undoubtably drives the observed jets 
and Herbig-Haro outflows from HD 163296. 

This scenario also requires a stellar magnetic field. For 
steady-state accretion, the disk is disrupted by the stel¬ 
lar magnetic field at a radius, Rm, roughly the corota¬ 
tion radius where the Keplerian angular velocity in the 
disk equals the angular velocity at the stellar surface (e.g., 
Shu et al. 1994). For HD 163296, this corresponds to 
Rm =10.67?* (adopting usim ~120 km s _1 as tabulated by 
van den Ancker et al. 1998). This radius is where the 
torque exerted by the magnetic field on the disk approxi¬ 
mately equals the viscous torque in the disk and is given by 
the Alfven radius, r*t, for a dipole field. The Alfven radius 
can be expressed in terms of the magnetic field strength at 
the stellar surface, the stellar mass, radius, and accretion 
rate (e.g., Frank, King, & Raine 1985). For HD 163296, 

r*i=7.2xl0 9 M^ 2 // 7 i?*^. Assuming ta=Rm results in a 

value of .B^IOM/ kG for the surface field strength for 
HD 163296. 

For the minimal accretion rate of Af 7 ~6xl0“ 5 , the re¬ 
quired magnetic field strength is only ~80 G. Higher mag¬ 
netic fields, of order 500 G, have been invoked to explain 
confinement of plasma during strong flares observed from 
the Herbig Ae star V892 (Giardino et al. 2004) and the 
Herbig Be star MCW 297 (Hamaguchi et al. 2000). Very 
recently, Hubrig, Scholler, & Yudin (2004) have presented 
evidence for (longitudinal) fields of several hundred Gauss 
on Herbig Ae stars — values consistent with surface field 
strengths of up to 2 kG measured from many classical 
T Tauri stars (e.g., Guenther et al. 1999; Johns-Krull, 
Valenti, & Koresko 1999). Theoretically, these are plau¬ 
sible relic field strengths for a young star (see references 
in the Introduction) and imply accretion rates of order 
AL 7 ~10 “ 3 or higher. Combining this with the simplistic 
arguments based on 1-D shocks given above suggests an 
accretion rate of Af~10 “ 9±1 M 0 yr -1 , tuned such that the 
great majority of the X-ray emission is reprocessed within 
the photosphere and such that only of order 1 % or less is 
actually observed at X-ray wavelengths. This rate is com¬ 
parable to those reported in the literature. Gullbring et 
al. (1998) report rates as low as a few 10 “ 10 M 0 yr -1 for 
a number of T Tauri stars while Calvet et al. (2004) find 
rates of a few 10 “ 8 M 0 yr -1 for young 1-3 M 0 T Tauri 
(pre-Herbig Ae) stars. 

The accretion rate and magnetic field strength of 
HD 163296 can be better determined when the contri¬ 
butions from the observed mass outflow are properly ac- 
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counted for in the overall energy budget. Until then we 
conclude that the X-ray emission from HD 163296 is most 
consistent with an optically-thin shocked accretion flow 
confined by a stellar (dipole) magnetic field of unknown 
origin. The X-ray data do not require a a late-type com¬ 
panion be present, and are in fact inconsistent with the 
hotter coronal temperatures (10 7 K) that typically char¬ 
acterize active late-type stellar coronae. A high accre¬ 
tion rate and strong magnetic field are also not required, 
though would ease difficulties that are otherwise present 
in explaining the accretion-driven X-rays using simplistic 
accretion shock models. 

HD 163296 is remarkable in being in a class of a very 
few pre-Main Sequence stars whose observed X-ray emis¬ 
sion has such a cool characteristic temperature that it 
can be entirely explained in terms of accretion. Accreting 
T Tauri stars generally exhibit X-ray emission character¬ 
istic of plasma at tempertures of 10' K or higher, similar 
to very active late-type stars. There is to our knowledge 
only one other PMS object whose X-ray emission is also 
accretion-dominated: the nearby T Tauri star TW Hya 
(e.g. Kastner et al. 2002). Other phenomena such as giant 
outbursts have also been observed (Kastner et al. 2004) 
that signal the onset of mass accretion at rates much higher 
than required to explain the X-ray behavior reported here 
for HD 163296. 

4.2. X-ray Knot 

Detection of X-ray emission associated with Herbig- 
Haro objects is not unprecedented. These shock-heated 
objects are formed when supersonic outflow collides with 
ambient material or previous ejecta (e.g., Schwartz 1983). 
Pravdo et al. (2001) report X-ray emission from HH2 
in Orion at the level of 5 x 10 29 ergs s” 1 and Favata et 
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al. (2002) estimate a luminosity of ~3 x 10 29 ergs s -1 
from HH knots in the proto-stellar jet of L1551 IRS5. 
The luminosity deduced here for the isolated knot in 
the SW jet is far less than these values yet is not 
unreasonable for an adiabatic shock. From the ana¬ 
lytic model of Raga et al. (2002), the luminosity is 
L-1.8xl0 29 (n/100) 2 (r/10 16 ) 3 (v/100) ergss" 1 for a shock 
velocity v km s -1 , a pre-shock density n cm -3 , and a bow 
shock length scale r cm. For a shock velocity equal to 
that of the jet (~350 km/s), a typical pre-shock density 
(~100 cm -3 ), and a bow shock length scale of order the 
optical size of the knot (~10 15 cm), the free-free luminos¬ 
ity is ~6 x 10 26 ergs s _1 . This is comparable to the value of 
~3xl0 27 ergs s” 1 estimated here for the emission feature 
along the HD 163296 jet axis. 
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